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Abstract 

The application of optical lattices allows a tuning of the geometry of Bose-Einstein condensates 
to effectively reduced dimensions. In the context of solid state physics the consideration of the 
__i . low-dimensional Frohlich polaron results in an extension of the polaronic strong coupling regime. 

^ I With this motivation we apply the Jensen-Feynman variational principle to calculate the ground 
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state properties of the polaron consisting of an impurity in a Bose-Einstein condensate in reduced 
dimensions. Also the response of this system to Bragg scattering is calculated. We show that 



'^ I reducing the dimension leads to a larger amplitude of the polaronic features and is expected to 

O ' facilitate the experimental observation of polaronic properties. In optical lattices not only Feshbach 

resonances but also confinement-induced resonances can be used to tune the polaronic coupling 
^ , strength. This opens up the possibility to experimentally reveal the intermediate and strong 

o 

CN I polaronic coupling regimes and resolve outstanding theoretical questions regarding polaron theory. 
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I. INTRODUCTION 



In recent years ultracold atomic systems have revealed themselves as quantum simulators 
for many-body theories [l|. Especially their high degree of tunability makes them attractive 
for this purpose. An example of a system that can be simulated in this way is the Frohlich 
polaron which is well-known from solid state physics where it is used to describe charge 

[i 

carriers in a polar solid (see for example Ref. [.2| for an extended overview). In the context 
of ultracold gases the system of impurities embedded in a Bose-Einstein condensation can 
be mapped onto the Frohlich polaron Hamiltonian |3|, |j] . In this case the role of the charge 
carriers is played by the impurities and the lattice vibrations are replaced by the Bogoliubov 



excitations. Recently this system has gained much interest both theoretically |5l-ll4| and 
experimentally 15l-ll9|. 

For the present work we focus on a single Frohlich polaron for which the Hamiltonian can 
not be analytically diagonalized and one has to rely on approximation methods. The most 
advanced theory for the ground state properties is the Jensen- Feynman variational princi- 
ple 20| which can be extended through the Feynman- Hellwarth-Iddings-Platzman (FHIP) 



approximation for the response properties 2l|, |22|] . The optical absorption of the Frohlich 



solid state polaron was later also obtained through a diagrammatic Monte Carlo calculation 
and a comparison with the FHIP approximation showed a good agreement at weak and 
intermediate polaronic coupling but in the strong coupling regime deviations were revealed 
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24| . Since there is no known material that exhibits the strong coupling behavior only 



the weak and intermediate coupling regime could be experimentally probed which resulted 
in a good agreement with the theory 25|, |26|]. A better understanding of the strong coupling 
regime could also shed light on the possible role of polarons and bipolarons in unconventional 



pairing mechanisms for high-temperature superconductivity [27|, |28|. Recently it was shown 
that for an impurity in a condensate the use of a Feshbach resonance allows an external 
tuning of the polaronic coupling parameter which makes it a promising system to probe the 
strong polaronic coupling regime for the first time [10|. Recent experiments have shown the 
feasibility of using Feshbach resonances for the tuning of interparticle interactions between 



29141 



different species |29H31|. 

Since the impurities are considered as not charged it is not possible to conduct optical 
absorption measurements to reveal the polaronic excitation structure as is possible for the 



Frohlich solid state polaron. It was shown in 13| that Bragg spectroscopy is suited to 
experimentally probe the polaronic excitation structure of an impurity in a condensate. 
Bragg scattering is a well established experimental technique in the context of ultracold 



gases (see for example Refs. [32 



4 l33|). The setup consists of two laser beams with different 
frequencies Ui and U2 and different momenta ki and k2 that are radiated on the impurity. 
The impurity can then absorb a photon from laser 1 and emit it to laser 2 during which 
process it has gained an energy u = ui — U2 and a momentum k = ki — k2. The response 

— * — * 

is reflected in the number of impurities that have gained a momentum k as a function of k 
and u. This number is proportional to the imaginary part of the density response function 
x(oo,k) [34]: 
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with pj: the density operator of the impurity. 

Another powerful tool in the context of ultracold gases is the ai 
which can be employed to modify the geometry of the system [35| . This allows to confine 
the system in one or two directions such that the confinement length is much smaller than 
all other typical length scales which results in an effectively low dimensional system. For 
these systems the interparticle interactions can be described through a contact pseudopo- 
tential with an amplitude that is a function of the three-dimensional scattering length and 
the confinement length. This permits to experimentally tune the interactions between the 
particles by varying the strength of the confinement which results in a resonant behav- 
ior. These confinement-induced resonances have been studied both theoretically 36|-|40| and 
experimentally 4l|-|45|. 

In the present work we adapt the calculations of the ground state and response properties 
of the polaronic system consisting of an impurity in a condensate to the case of reduced 
dimensions. This was done for the Frohlich solid state polaron in Refs. 46l448| which led 



to the polaronic scaling relations. These are applicable for polaronic systems of which the 
interaction amplitude V^ (see later) is a homogeneous function. Unfortunately this is not the 
case for the polaron consisting of an impurity in a Bose-Einstein condensate. We start by 
showing that also in lower dimensions the Hamiltonian of an impurity in a condensate can 
be mapped onto the Frohlich polaron Hamiltonian. Then the Jensen-Feynman variational 
principle is applied to calculate an upper bound for the free energy and an estimation of the 



effective mass and the radius of the polaron as was done in Ref . [10| for the three-dimensional 
case. Subsequently the treatment of Ref. [l3y for the response to Bragg spectroscopy in 3 
dimensions is adapted to reduced dimensions. All results are applied to the specific system 
of a lithium-6 impurity in a sodium condensate. 

II. IMPURITY IN A CONDENSATE IN d DIMENSIONS 

The Hamiltonian of an impurity in an interacting bosonic gas is given by: 

k k,k',q k,q 

The first term in this expression represent the kinetic energy of the impurity with p {r) 
the momentum (position) operator of the impurity with mass mj. The second term in the 
right-hand side of ([2]) describes the kinetic energy of the bosons with creation (annihilation) 

operators \at> ({%}) and energy E^ = ^ fi where /i is the chemical potential of the 

bosons and ms their mass. The last two terms represent the interaction energy with Vbb (<?) 
the Fourier transform of the boson-boson interaction potential and Vjb {(f) of the impurity- 
boson interaction potential. All vectors in expression ([2]) are considered as (i-dimensional. 

In Refs. J49| and 50| it is shown that in one and two dimensions, respectively, at tem- 
peratures well below a critical temperature T^ a trapped weakly interacting Bose gas is 
characterized by the presence of a true condensate while just below Tc this is a quasiconden- 
sate. A quasicondensate exhibits phase fluctuations with a radius R^ that is smaller than 



the size of the system but greatly exceeds the coherence length ^ [49|, |50| . Since the radius of 
the polaron Rpoi is typically of the order ^ (see later) we have Rpoi <C -R^ which shows that 
the polaronic features are also present in a quasicondensate. In the following we no longer 
make the distinction and use the name condensate for both situations. The presence of a 
condensate can be expressed through the Bogoliubov shift which (within the local density 
aproximation) transforms the Hamiltonian ([2]) into [10|: 

H = EGP + giBNo + Hpa, (3) 

where use was made of contact interactions, i.e. Vbb {(f) = Qbe and Vib [q] = Qib- In order 
to have a stable condensate the boson-boson interaction should be repulsive, i.e. qbb > 0. 
The sign of the impurity-boson interaction strength gjB is in priciple arbitrary, however 



for the Bogoliubov approximation to be valid the depletion of the condensate around the 



impurity must remain smaller 



;he condensate density which means the formalism is not 



;han 
valid for a large negative gjB M 1^ • The first term in the right-hand side of (^ represents 
the Gross Pitaevskii energy Ecp of the condensate and the second term gives the interaction 
of the impurity with the condensate (with A'"o the number of condensed bosons in a unit 
volume). The third term is the polaron Hamiltonian which describes the interaction between 
the impurity and the Bogoliubov excitations: 

Hpoi = ^ + 5Z ^fc"|% + Yl ^kPi (^) (% + "Lfe) ' (4) 

where \ot\ ({%}) are the creation (annihilation) operators of the Bogoliubov excitations 
with dispersion: 



with ^ the healing length: ^ = \/ 2m n • '^^^ interaction amplitude V^ is given by: 

'-v'iVo.,.!^)'". (6) 

III. POLARONIC GROUND STATE PROPERTIES IN d DIMENSIONS 

In this section we summarize the main results from standard polaron theory regarding 
the ground state properties with emphasis on the dependency on the dimension (see for 



example Ref. 52| for more details) and we apply this to the polaronic system consisting of 



an impurity in a condensate. 



A. Jensen-Feynman variational principle 



The most accurate available description of the ground state properties of a polaron is 



based on the Jensen-Feynman inequality which states [20|, l53| : 



J'<^o + ^{S-So),^, (7) 

with J^ the free energy of the system, J^q the free energy of a trial system, /3 = (ksT)' 
the inverse temperature with ks the Boltzmann constant, S the action of the system and 



So the action of the trial system. It was suggested by Feynman to consider the particle 
harmonically coupled to a mass M with a coupling constant MW'^ for the trial system 20|. 
This leads to the following expression for the Jensen- Feynman inequality ([7]) 46l |: 
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with d the dimension, V the volume, fi = W \/\ + M/mi and Q \k^u\ the Green function 
of the Bogoliubov excitations: 

cosh [uj^ {u - hP/2)] 



Q {k,u 
and M.M,n {k,u] the memory function: 

he 



sinh [h(3uj^/2] 
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The parameters Q and M are then determined variationally by minimizing the expression 
(jH]). The present treatment also allows an estimation of the radius of the polaron as the root 
mean square of the reduced coordinate r of the model system 54l |: 



, 2\ , h nil + M , f(3hQ 
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In 20| Feynman also presented a calculation of the polaronic effective mass m* at zero 
temperature: 
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2 (m + M) fi 

As far as we know there exists no generalization of equation ( IT2l) to finite temperatures but 
as a first estimation we use (fT2|) with the temperature dependent variational parameters M 
and VL. 



B. Polaron consisting of an impurity in a condensate 



Here we introduce the Bogoliubov spectrum ([5]) and the interaction amphtude IQ which 
are specific for the polaronic system consisting of an impurity in a condensate. This allows 
us to write the Jensen-Feynman inequality ([H]) as (we also use polaronic units, i.e. h = ^ = 
mi = 1): 
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where we introduced the dimensionless coupling parameter a'^'^'^ as follows: 
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with r (x) the gamma function. The prefactor was chosen to be in agreement with the 
definition for a^^-* in Ref. [10| ■ Note that the coupling parameter depends on the impurity- 
boson interaction amplitude qjb and also on the boson-boson interaction amplitude gsB 
through the healing length C,- As mentioned in the introduction these interaction amplitudes 
and thus also the coupling parameter can be externally tuned through a Feshbach resonance 
or in reduced dimensions also with a confinement induced resonance. 

For d = 2 the fc-integral in flT4l) contains an ultraviolet divergence. This is also the case 

n 

in 3 dimensions and it was shown in Ref. [10] that this is solved by applying the Lippmann- 
Schwinger equation up to second order for the interaction amplitude in the second term of 
the Hamiltonian (131). This results in a renormalization factor that is incorporated through 

n 

the following substitution [10|: 
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with T (E) the scattering T-matrix. In 2 dimensions the limit i? — )■ in ( lT6l) results in an 
infrared divergence. The second term in (lT6ll can be written as: 

1 a(2) fi^ mB + mi f°° k 
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which hfts the ultraviolet divergence in (fT^ . For numerical considerations a cutoff Kc is 
introduced for the fc- integral which enables us to calculate the integral in flT7|) : 

a^^) h? rriB + mi f^" k ,, a^^^ ff niB + mr ( ^^ — E 

:dk = -r m 



27r r/i/^^ r/i^ Jg /c^ — 2mr^'^E/fi^ Att rrij^'^ tjib \ E 

^o^ffm, + „.,JrfKl_Y (18) 

A-K mii^ TUB \2mrEj ^ ^ 

where in the second step we used the fact that the energy related to the cutoff is much larger 
than the typical energy of the scattering event E. Equation ( TT8l) shows that the chosen value 
of E is not important since it only results in an energy shift and therefore has no influence 
on the physical properties of the system. 

C. Results 

We apply the presented treatment to the system of a lithium-6 impurity in a sodium 
condensate {niB/mi = 3.82207). All results are presented in polaronic units, i.e. h = ^ = 
mi = 1. 

In figure [U the results for the polaronic ground state properties in 2 dimensions as a 
function of the coupling parameter a*^^-* are presented. In (a) the radius of the polaron is 
shown and in (b) the effective mass at different temperatures. The observed behavior is 
analogous to the three-dimensional case (see Ref. 10|]) and suggests that for growing a^'^^ 



the self-induced potential becomes stronger leading to a bound state at high enough a^'^K 
However, as compared to the three-dimensional case, the transition is much smoother with 
a transition region between a*-^-* Ri 1 and a^^' ~ 3. This behavior is in agreement with the 

□ n 

mean- field results of Refs. [7|, |55[, where also a smooth transition to the self-trapped state 
was found for d = 2. For the cutoff Kc we used the inverse of the Van der Waals radius 
for sodium which results in Kc = 200. To check whether this cutoff is large enough the 
variational parameter M is plotted in the inset of figure ([I])(b) for different values of Kc 
which reveals already a reasonable convergence at ii"c ~ 5. 

In figure |2] the results for the 1-dimensional case are presented. In (a) the radius of the 
polaron is plotted and (b) shows the effective mass at different temperatures. For growing 
a*-^-* the characteristics of the appearance of a bound state in the self- induced potential are 
again observed. The characteristics of the weak coupling regime are however not present 

8 




FIG. 1: (Color online) The groundstate properties of the polaron consisting of a lithiuni-6 impurity 
in a sodium condensate in 2 dimensions. In (a) the radius of the polaron (jlip is presented and 
in (b) the effective mass (J12p as a function of the polaronic coupling parameter a^^' at different 
temperatures (/3 = {ksT)^ ) and with a cutoff Kc = 200. The inset shows the variational param- 
eter M for different values of the cutoff Kc ai j3 = 50. All results are presented in polaronic units 
(h = mi=i=l). 

and the transition region is between a^^^ = and a^^^ ~ 1. This is asrain in agreement with 

nn 

the mean-field resuhs of Refs. p, |55[ for rf = 1. 

IV. RESPONSE TO BRAGG SCATTERING IN d DIMENSIONS 



The response of a system to Bragg spectroscopy is proportional to the imaginary part of 
the density response function ([1]). In Ref. 13| it was shown that the use of the Feynman- 
Hellwarth-Iddings-Platzman approximation (as introduced in Ref. 2l|] for a calculation of 
the impedance of the Frohlich solid state polaron and generalized in Ref. [22| for the optical 




FIG. 2: (Color online) The polaronic ground state properties of a litliium-6 impurity in a sodium 
condensate in 1 dimension. The radius (a) and the effective mass (b) are presented as a function 
of the polaronic coupling parameter a^'^' at different temperatures (/3 = {ksT)" ). All results are 
presented in polaronic units {h = mi = ^ = 1). 



absorption) leads to the following expression for the Bragg response: 
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For numerical calculations the representation for the self energy as derived in appendix |A] is 
used. 
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A. Sum rule 



As was first noted in [56| for tlie Frolilicli polaron and generalized in 13| for an impurity 
in a condensate the f-sum rule can be written as: 
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with e a small number such that the Drude peak (see later) is not included in the integral 



and: 
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In the limits /3 — )■ oo and fc — )■ the function (!24l) is related to the Feynman effective mass 
dnD 56|: 

(25) 

rom the 



m* = mj 1 — lim R (a, 0) 
This relation provides a powerful experimental tool to determine the effective mass 



optical response which was recently applied for the Frohlich solid state polaron 57 



B. Self energy for an impurity in a condensate 



Introducing the interaction amplitude ([6]) and the coupling parameters (1151) in expressions 
IA9|) and ( ]A12|) for the imaginary and real part of the self energy results in (using polaronic 
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See appendix|X]for the definition of tlie different functions. Tliese expressions are suited for 
numerical calculations of the Bragg response. 

C. Weak coupling limit 



At weak polaronic coupling the Bragg response flT9l) to lowest order in a is given by (in 
polaronic units): 

Im \x^ (co, k)] = -— Im [s^ fw, k)] . (28) 
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In the weak coupling limit the variational parameter M tends to zero and the imaginary 
part of the self energy fl26|) becomes: 
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These expressions coincide with the weak coupling result obtained in the framework of 
Gurevich, Lang and Firsov 59|. 



D. Results 

We present the Bragg response for a lithium-6 impurity in a sodium condensate 
{niB/mi = 3.82207). All results are presented in polaronic units, i.e. h = C, = frij = 1. 

In figure |3] the Bragg response (iT9l) is presented for different temperatures and for a 
momentum exchange A; = 1 in 1 and 2 dimensions at weak polaronic coupling (a^^^ = 0.01 
and a^^-* = 0.1). In both cases we observe the Drude peak centered at a; = and a peak 
corresponding to the emission of Bogoliubov excitations. This is qualitatively the same 



behavior as in the three-dimensional case 



13| . quantitatively we observe that the amplitude 



of the Bogoliubov emission peak increases as the dimension is reduced. The Drude peak is 



a well- 
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cnown feature in the response spectra of the Frohlich polaron (see for example Refs. 



60, 



6l| ) and is a consequence of the incoherent scattering of the polaron with thermal 



Bogoliubov excitations. The width of the Drude peak scales with the scattering rate for 

absorption of a Bogoliubov excitation which is proportional to the number of thermally 

i 
excited Bogoliubov excitations [62] . This explains the temperature dependence of the width 

of the Drude peak in figure |3l 

In ID another sharp peak is observed in figure [3] at u; = 0;^ (with Uk the Bogoliubov 
dispersion ([5]) and k the exchanged momentum) which broadens as the temperature is in- 
creased and dominates the Bogoliubov emission peak at relatively high temperatures. This 
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FIG. 3: (Color online) The Bragg response (I19p at weak polaronic coupling, momentum exchange 
k = 1 and for different temperatures (/3 = {ksT)^ ) in ID (a) and 2D (c). In both cases a peak 
corresponding to the emission of Bogoliubov excitations is observed together with the anomalous 
Drude peak at cj = 0. In ID another sharp peak is present at uj = uJk, with cok the Bogoliubov 
dispersion ([5|). In (b) we have zoomed in on this sharp peak in ID. All quantities are in polaronic 
units {h = mi = ^ = 1). 



extra peak in ID is associated with the weak couphng regime since at intermediate couphng 
the sharp structure disappears and the peak merges with the Bogohubov emission peak. 
The location indicates that it corresponds to the process where both the exchanged energy 
huj and momentum k are transferred to a Bogoliubov excitation. Whether this extra peak 
is experimentally observable is questionable since it is only visible at relatively high temper- 
atures where in reduced dimensions thermal phase fluctuations can become important and 
destroy the polaronic features. 

Figure S] presents the Bragg response for different momenta exchange at a temperature 
/3 = 100 (where the sharp peak at the Bogoliubov dispersion in ID is too narrow to perceive). 
The insets show the location of the maximum of the Bogoliubov emission peak as a function 
of the exchanged momentum together with a least square fit to the Bogoliubov spectrum 
(|5]) which results in a good agreement. The optimal fitting parameter is determined as 
niB = 4.3159 (4.2216) in ID (2D). 

In figures [5] and [6] we have zoomed in on the tail of the Bogoliubov emission peak for dif- 
ferent values of the coupling parameter in 1 and 2 dimensions, respectively. At larger values 
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FIG. 4: (Color online) The Bragg response at weak polaronic coupling {a^ ' = 0.1) for different 
exchanged momenta k in 1 and 2 dimensions and at a temperature /3 = 100. The inset shows the 
location of the maximum of the peak as a function of the exchanged momentum (markers) together 
with a least square fit to the Bogoliubov spectrum ([5]) (full line), this results in ms = 4.3159 
(4.2216) for the fitting parameter in ID (2D). Everything is in polaronic units {h = rrij = ^ = 1). 



for the polaronic coupling parameter a^'^'^ the emergence of a secondary peak is observed. 
This behavior is also observed in the optical absorption of the Frohlich solid state polaron 
where the secondary peak corresponds to a transition to the Relaxed Excited State accom- 
panied by the emission of phonons [63]. The Relaxed Excited State denotes an excitation of 
the impurity in the relaxed self-induced potential where relaxed means that the self-induced 
potential is adapted to the excited state wave function of the impurity. In the inset the lo- 
cation of this secondary peak is plotted as a function of the exchanged momentum together 
with a least square fit to the following quadratically spectrum: 



00 {k) = u + 



2m 



(31) 



which shows a good agreement at small k. This suggests that the state corresponding to the 
secondary peak is characterized by a transition frequency u and an effective mass m (this 
was also observed for the 3-dimensional case in Ref. 13|). 

Finally we have checked whether the spectra satisfy the sum rule ( 123|) . We calculated the 
sum of the two terms on the left hand side of expression fl23l) which is presented in table [T] 
for d = 1 and in table [ITl for (i = 2 at /3 = 100 and at different values for a and k. These 
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FIG. 5: (Color online) Here we zoomed in on the tail of the Bogoliubov emission peak for momentum 
exchange k = 1 and temperature /3 = 100 in 1 dimension. It is clear that at larger values for a^^' 
a secondary peak emerges. The inset shows the location of the maximum of this secondary peak 
at a^^' = 3 as a function of the exchanged momentum (markers) together with a least square fit 
to a quadratic spectrum ([3T|) (solid line), this results in a; = 1.6386 and m = 2.0107 for the fitting 
parameters. Everything is in polaronic units {h = mj = .^ = 1). 




FIG. 6: (Color online) Here we zoomed in on the tail of the Bogoliubov emission peak for momentum 
exchange k = 1 and temperature f3 = 100 in 2 dimensions. As in the one-dimensional case a 
secondary peak emerges at larger values for a^'^'. The inset shows the location of the maximum 
of this secondary peak at a^'^' = 4 as a function of the exchanged momentum (markers) together 
with a least square fit to a quadratic spectrum (|3ip (solid line), this results in a; = 2.0601 and 
m = 2.0755 for the fitting parameters. Everything is in polaronic units {h = mj = (, = 1). 
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k = 1 

k = 3 



a 



(1) = 0.1 qW = 3 



1.5440 
1.5544 



1.5547 
1.5743 



TABLE I: Here we show the sum of the two terms at the left hand side of the f-sum rule (j23p in 1 
dimension at /3 = 100 and at different values for a^^' and k. 



k 
k 



a 



(2) 



1 a 



(2) 



1.5678 1.5734 
1.5669 1.5800 



TABLE II: Here we show the sum of the two terms at the left hand side of the f-sum rule (j23p in 
2 dimensions at /? = 100 and at different values for a^"^' and k. 

values should be compared to 7r/2 = 1.5708 which results in a fair agreement with small 
deviations which are to be expected since numerically we had to introduce a cutoff for the 
w-integral in ( 123|) and the choice of the parameter e in ( l23l) is somewhat arbitrary resulting 
in a double counting of part of the weight of the Drude peak. 

V. DISCUSSION AND CONCLUSIONS 

We have applied the calculations for the polaronic ground state properties of an impurity 
in a Bose-Einstein condensate and the response of this system to Bragg spectroscopy to 
reduced dimensions. For this purpose we introduced a polaronic coupling parameter a^'^^ ( !T5|) 
which depends on the dimension. For growing a^"^^ the ground state properties suggest that 
the self-induced potential accomodates a bound state. As compared to the three-dimensional 
case the transition to the self-trapped state is much smoother in reduced dimension and for 
d = 1 the characteristics of the weak coupling regime are absent. 

The Bragg response of the system revealed a peak corresponding to the emission of 
Bogoliubov excitations, the Drude peak and the emergence of a secondary peak in the 
strong coupling regime. The amplitude of these polaronic features grows when we go to 
reduced dimensions. This is important since this indicates that going to reduced dimensions 
can facilitate an experimental detection of polaronic features. In ID another sharp peak is 
observed at weak polaronic coupling that corresponds to the full transition of the exchanged 
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energy and momentum to a Bogoliubov excitation. 

Another advantage of considering reduced dimensions is the possibihty of using 
confinement-induced resonances which permits a tuning of the polaronic couphng param- 
eter. These resuhs show that considering an impurity in a Bose-Einstein condensate in 
reduced dimensions is a very promising candidate to experimentally probe the polaronic 
strong coupling regime for the first time. 
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Appendix A: Other representation for the self energy 



Here we rewrite the self energy fl20|) to a form which is more suited for numerical calcu- 
lations. The presented derivation is based on the approach for the optical absorption of the 



Frohlich solid state polaron as proposed in Refs. [22|, |6J]. We start by rewriting D (t) ( 122|1 
as: 

t^ . h 

^ ' ~ 2/3 (m + M) ~ *2(m + M) 

which allows us to write: 



with: 



a^ {13) = — -— coth — — ; 

B{(3,n,n') = ^^^{a^l + nmy[a'n{n)f; (A3) 
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and a = a (oo). If we now use (lA2p in the expression for the self energy (l20l) we get: 

2 

2(n+n') 



T,[uj,k 



9 ^ 1 /i;.o 



miNh 



n,n'=0 qj^O 



A;2 



/;; + g 



S(/3,n,n')e-"'(^K^>9) 



X 



f rft(l-e''"*)lm|[l + nK-)] 



(fc+.r)^ 



ri^+i^+t 



(k+?f_^2 



g 2,9(m+M)'' ^'■^' " _|_ ^ (^^-.j g 2/3(m+M) 



+iA^t 



with: 



k + q) h 

A± = ±ujff+ ^- ^ + fi (n - n') . 



(A4) 



(A5) 



"^ ' 2(m + M) 

We now spht the self energy in an imaginary and a real part. Taking the imaginary part of 
H) results in: 



Im 



T, (uj,k 



2 v^ V-,.. 2 ^^"^ 



2{n+n') 



miNh 



E Ei^^S^ (^^+^T^ ^^^'"'"^"^^''^^'"'^ 



n,n'=0 (fT^O 



X Im [ / dt sin (wt) < [1 + n (a;,-)] e 2/3(-lAf)*'+*^+* 



(A6) 



( k-\-Q] 

The time integration can now be done (with C^ = 2g(m+Af ) ) • 



(it sin (wt) < [1 + n (wg-)] e 



-C^i^+.A+t^^(^^^)g-CV+.A-t 



1 f 


L/¥ K-)^ 


^A+ + cu\ .AF K-") 


V 2C ) 
\ 


2iC \ 


[1 + n {Ug-)\ 


^g 4C^ + iF 

T .^fA-+uj\ 


y 2C J 2 


+n (wg-) 


-^— e 4c^ 
2 


\ 2C J 


2 \ 2C J 


)• '' 



Where we introduced the Dawson integral F (x) : 



F(x) = e-" / e^ dy 
Jo 

= -y/ne'"^ erfi(x) 



(A8) 



where erfi(x) is the imaginary error function: erfi(x) = —ieii{ix), with erf (x) the error 
function. This finally results in the following expression for the imaginary part of the self 
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energy: 



Im 



Tj {u,k 



y^27r f3 (m + M) 



X e 



2mTNh 



-a^mk+q) 



EEi*-. 



k.q 



2(n+n')-l 



A; + g 5(/3,n,n') 



n,n/=0 <f^0 



[1 + ^K-)] 



A;2 



j3{m + M)(A+ +uj\ ^(m + M)(A+ -uiY 



g 2(ft + g)^ 



g 2(k + qf 



+n {UJ^) 



I3{m+M){A +ui\ /3(m+Af)f A^-ij) 

g 2(fe+9)^ _ g 



2{k+qf 



(A9) 



For the real part of flA4p we have: 



Re 



Tiioj.k 



EEi*-. 



rajNh ^-~ 

n,n'=0 g^O 






fc2 



X Im ( / dt [1 - cos (wt)] <j [1 + n (wg-)] e 2/*("'+Af) 



V+iA't 



t'^+iA+t 



+n(Uff)e a/JC^+w) 



(AlO) 



The time-integration is in this case: 



dt [1 - cos (cot)] |[1 + n (u^)] e-^'*'+^^+* + n (u^) g-f^'t'+*^-t| 



^<([l + nK-)] 






V .g-^j^_,ir(___ 

A' +L0 



+ n (ojg-) 



^e ^c^ + 2iF 



2C 



\A-+u 



e Jc^ — ii^ 



2C 



2 V 2^ 



(All) 
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This results in the following expression for the real part of the self energy: 

k-q] /^ x2{n+n')-l 



Re 



T. \uj,k 



v/2/3 (m + M) 
miNh 



T. j:k' 



n,n'=0 q=iO 

^/2/3 (m + M) (A+ + io) 



k^ 



k + q 



B{/3,n,n') 



^2/3 (m + M)A^ 



+ n (u^) { 2F 



k + q 

^2/3 (m + M)A 



k + q 

y/2/3 (m + M) (A+ - w) 



/c + (? 



-F 



/c + g 



v/2/3 (m + M) (A- + w) 



/c + g 



v/2/3 (m + M) (A- - to) 



k + q 



(A12) 
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